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Abstract. The paper concerns lattice triangulations, i.e., triangulations of the integer 
points in a polygon in R 2 whose vertices are also integer points. Lattice triangulations 
have been studied extensively both as geometric objects in their own right and by virtue 
of applications in algebraic geometry. Our focus is on random triangulations in which 
a triangulation a has weight A' CT ', where A is a positive real parameter and \a\ is the 
total length of the edges in a. Empirically, this model exhibits a "phase transition" 
at A = 1 (corresponding to the uniform distribution): for A < 1 distant edges behave 
essentially independently, while for A > 1 very large regions of aligned edges appear. We 
substantiate this picture as follows. For A < 1 sufficiently small, we show that correlations 
between edges decay exponentially with distance (suitably defined), and also that the 
Glauber dynamics (a local Markov chain based on flipping edges) is rapidly mixing (in 
time polynomial in the number of edges in the triangulation). This dynamics has been 
proposed by several authors as an algorithm for generating random triangulations. By 
contrast, for A > 1 we show that the mixing time is exponential. These are apparently 
the first rigorous quantitative results on the structure and dynamics of random lattice 
triangulations. 



1. Introduction 

1.1. Background. Let A^ n = {0, 1, . . . , rei} x{0, 1, . . . , re} denote the set of integer points 
in an m x n rectangle in M?. This paper is concerned with (full) triangulations of n , i.e., 
triangulations that use all the points. Any such triangulation partitions the rectangle into 
2mn triangles each of which is unimodular (has area 1/2). See Figure [T] for an example of 
a triangulation with m = 5, n = 7. Most of our results extend to the case of triangulations 
of the integer points in an arbitrary (not necessarily convex) lattice polygon. 

Lattice triangulations are fundamental discrete geometric objects that have a rich and 
beautiful structure (see, e.g., [3] for background). For example, a number of recent papers 
have developed ingenious combinatorial arguments to estimate the asymptotic number 
of triangulations [TJ JSJ \T3[ [Tfi] , as well as to explore their connectivity properties under 
natural local moves [8] . 

Lattice triangulations have also received much attention in algebraic geometry, through 
connections with plane algebraic curves and Hilbert's 16th problem [21j . the theory of 
discriminants [6] , and toric varieties [2] . In several of these contexts one is chiefly interested 
in properties of "typical" triangulations, such as whether they are regular (in the sense 
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(a) A = 1 (b) A = 1.1 (c) A = 0.9 



Figure 2. Random triangulations of a 50 x 50 lattice 



of being representable by a "nice" lifting function [2Q|, I24j). whether they contain "long, 
thin" triangles, etc. This leads one to investigate random triangulations drawn from a 
uniform distribution O [22] . 

A natural generalization that can yield useful insights here is to introduce weights: 
specifically, we consider the distribution in which a triangulation a has weight A' "', where 
\a\ is the total t\ length of the edges in a and A > is a real parameter. The case A = 1 
is the uniform distribution, while A < 1 (respectively, A > 1) favors triangulations with 
shorter (respectively, longer) edges. Figure [2] shows typical triangulations of a 50 x 50 
region for three values of A. This weighted version corresponds rather closely to a model 
that has recently been studied in statistical physics [TTJ. 

One striking feature of the pictures in Figure [2] is the tendency in the case A > 1 for 
edges to line up in macroscopically large regions of similar slope, while for A < 1 these 
regions disappear. The uniform case A = 1 appears to represent a "phase transition" 
between these two regimes, in which fairly large regions form but are unstable. One goal 
of this paper is to initiate a rigorous, quantitative study of this phenomenon. 

How does one generate a random triangulation such as those shown in Figure [2]? The 
only known feasible approach for large values of m, n is to simulate a local Markov chain, 
or Glauber dynamics, which converges to the desired probability distribution over trian- 
gulations. Fortunately there is a very natural dynamics here: pick a random edge of the 
triangulation and, if it is the diagonal of a convex quadrilateral (which is in fact always 
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a parallelogram), flip it to the opposite diagonal. The induced graph on the set of tri- 
angulations, in which two triangulations are adjacent iff they differ by one edge flip, is 
known as the "flip graph." The flip graph is well known to be connected [S] (indeed, quite 
a bit more is known about its geometry [5j[8]), so the Glauber dynamics converges to the 
uniform distribution over triangulations. Moreover, a standard "heat-bath" modification 
of the flip probabilities achieves the weighted distribution A' "' for any desired A. 

This leads to a fundamental question of both algorithmic and structural interest, posed 
by Emo Welzl and others |22| : what is the mixing time of the Glauber dynamics, i.e., 
the number of random flips until the dynamics is close to the stationary distribution? 
The pictures in Figure El together with the actual Glauber dynamics simulations used to 
generate them, suggest the following conjecture. When A < 1 there are no long-range 
correlations between edges, and we would expect the mixing time to be small; when 
A > 1, the large "frozen" regions of aligned edges take a very long time to break up, and 
consequently the mixing time should be large. We may summarize this in: 

Conjecture 1.1. (a) For any fixed A < 1, the mixing time of the Glauber dynamics is 
poly(m + n) (specifically, 0(mn(m + n))). (b) For any fixed A > 1, the mixing time is 
exponential in (m + n) (specifically, exp(J7(mn(m + n))) ). 

The "critical" case A = 1 appears to be rather delicate: while large regions tend to 
form, they tend not to be stable over long time scales, and the behavior of the mixing 
time is less clear. 

A second goal of this paper is to make the first rigorous progress in analyzing the mixing 
time of these dynamics, and in particular towards a proof of Conjecture 11.11 

Before describing our results, we relate random triangulations to classical spin systems in 
statistical physics and explain why they present a challenge to existing analysis techniques. 
It is well known that, in any triangulation of A^ n , the midpoints of all the edges are 
fixed, and are precisely the points of the half-integer lattice (with the original lattice 
points removed); see the example in Figure [TJ Thus we can view any triangulation as 
an assignment of "spin values" (edges) to each of these midpoints, subject of course to 
consistency constraints. One might hope, then, to capitalize on the vast literature on both 
the structure and dynamics of spin systems. However, there is a crucial difference here. 
In a classical spin system on a graph, all interactions are local; i.e., the distribution of 
the value of a given spin is determined by the spins of its neighbors in the graph. For 
triangulations, on the other hand, each edge belongs to two triangles and thus has four 
neighboring edges — but the midpoints of these edges depend on the triangulation and 
may lie very far from the midpoint of the edge itself. In other words, the geometry of 
a triangulation seems to have little to do with the Cartesian geometry of the lattice. It 
is this non-locality that makes the application of established techniques very challenging 
here. 

1.2. Contributions. Our first results confirm the empirical observation that, in the 
regime A < 1, correlations between edges decay rapidly with distance — a property of- 
ten referred to in the spin systems literature as "strong spatial mixing." Specifically, for 
any sufficiently small A < 1 , we show that if the configuration of the edge at midpoint x is 
fixed to a x , then the effect on the distribution of the edges at other midpoints z £ A, for 
an arbitrary subset A, decays exponentially with a natural distance d(A,a x ) between A 
and a x . (See Section U] for the definition of d(-, •).) 
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Theorem 1.2. There exists Ao > such that for all A < Ao the following holds. Let 
H denote the X-weighted distribution on triangulations o/A^ n , and [f^ the distribution 
conditional on the edge at x having configuration a x . Then the variation distance between 
[i and jjL ax at any set of midpoints A satisfies — /x^Ha ^ \A\ exp(— cd(A, o~ x )), where 
the constant c > depends only on A. 

We remark that this result holds in the presence of arbitrary fixed edges, or constraints, 
and therefore for triangulations of arbitrary lattice polygons. Handling constraints in- 
troduces some technical complications into our proofs, and necessitates in particular the 
understanding of minimum length, or "ground state" triangulations. (In the absence of 
constraints, ground state triangulations are trivial: all edges are either horizontal, vertical 
or unit diagonal.) Fortunately, as we show, ground state triangulations with constraints 
can be constructed greedily, by placing each edge independently in its minimum length 
configuration consistent with the constraints. We also prove another fundamental property 
in this regime, namely that the probability of an edge exceeding its ground state length 
by k decays exponentially with k. 

Theorem 11.21 is analogous to spatial mixing results in classical spin systems, with the 
additional twist of the distance function d whose definition is tailored to the geometry 
of triangulations. Our argument is reminiscent of classical "Peierls-type" arguments for 
the Ising model, which explains why we require A to be sufficiently small. We conjecture 
that the theorem holds for all A < 1, but handling values of A close to 1 will require more 
sophisticated methods. 

We then turn to the Glauber dynamics, again for the regime A < 1. Here we can show 
the following. 

Theorem 1.3. There exists \\ > such that for all A < Ai the mixing time of the 
Glauber dynamics on triangulations of n with parameter A ( and in the presence of 
arbitrary constraints) is 0(mn(m + n)). 

Theorem 1 1 . 31 goes some way towards proving part (a) of Conjecture II .1| though it does 
require A to be sufficiently small (rather than merely less than 1). Our proof here uses 
a path coupling argument with a suitably chosen exponential metric on triangulations, 
and exploits properties of the geometry of the flip graph which we also establish. We also 
prove that our bound on the mixing time is tight. In the special case of triangulations of 
a 1-dimensional region A® n , we show that the mixing time is 0(n 2 ) for all A < 1. In this 
case triangulations turn out to be isomorphic to "lattice paths" of length 2n starting and 
ending at with ±1 increments at each step. 

Our final main result concerns the Glauber dynamics in the complementary regime 
A > 1. Here we are able to prove a slightly weaker version of part (b) of Conjecture ll.il 

Theorem 1.4. For any A > 1, the mixing time of the Glauber dynamics on triangulations 
of n with parameter A is exp(0(m + n)). 

While this result confirms the conjectured exponential slowdown for all A > 1, we believe 
that the mixing time should be exponential in the maximum total edge length mn(m + n) 
rather than in the maximum length of a single edge (m + n) . Theorem 11.41 is proved by 
exhibiting an explicit bottleneck in the dynamics, specifically an initial configuration from 
which it takes a very long time to change the orientation of some long edge. We also prove 
a stronger lower bound exp(0(n 2 /m)) which holds whenever m -C yfn. 

1.3. Related work. Triangulations of general point sets are a large topic with numerous 
applications in mathematics and computer science, including combinatorics, optimization, 
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algebraic geometry, computational geometry and scientific computing. For excellent back- 
ground the reader is referred to the survey by Lee [10] and the recent book of DeLoera, 
Rambau and Santos [3J, which also discusses lattice triangulations in some depth. As men- 
tioned earlier, lattice triangulations have been studied both in their own right as geometric 
objects, and in several contexts in algebraic geometry [2j l6l 120 ] [21 ] 124]. 

Much of the work on lattice triangulations has focused on counting them. A sequence 
of beautiful combinatorial arguments |16[ [IJ EJ [13] have shown that the number of trian- 
gulations of k° m n is at most 0(6.86"™) [13] and at least ft(4.15 mn ) [8]. Although it is not 
the tightest upper bound known, we briefly mention the elegant result of Anclin [lj, who 
shows that if the edges of a triangulation of A^ n are added one by one starting at the 
top left and proceeding left-to-right and top-to-bottom (by midpoint), then the maximum 
number of choices for each edge is two. Since there are fewer than 3mn interior edges, this 
immediately yields an upper bound of 8 mn on the number of triangulations. (In contrast, 
the best known upper bound for the number of triangulations of a general set of n points 
in the plane is 43 n [18].) 

The flip graph on triangulations of n has also been studied in some depth, though we 
should stress that our work is apparently the first to handle constraints (and thus general 
lattice polygons). For example, the flip graph is known to have diameter 0{mn{m + n)) 
and to be isometrically embeddable into the hypercube [5]. Random walk on this graph 
has been proposed by several authors as an algorithm for generating random triangula- 
tions, and has been used heuristically in the formulation of conjectures regarding typical 
triangulations [HI [22] • However, nothing is known rigorously about its mixing time. (We 
note in passing that the analogous flip dynamics for triangulations of (the vertex set of) 
a convex polygon has been extensively analyzed [141 115] : however, that problem is much 
easier because the set of triangulations has a Catalan structure.) 

Random lattice triangulations with weights have appeared in slightly different form as 
a model in statistical physics [E]; there, triangulations are weighted according to the 
sum of squares of their vertex degrees. This is actually very close to our model, as high 
degrees are associated with long edges. The structural results reported in [T7] are based 
on simulations and are non-rigorous. Rigorous results have been obtained for a loosely 
related "topological glass" model by Eckmann and Younan [3]. 

The literature on structural properties and Glauber dynamics of lattice spin systems 
is too vast to summarize here. We refer the reader to the standard references [19] for 
structural properties such as spatial mixing, and |12] for mixing times of Glauber dynamics. 
As explained earlier, while triangulations may be viewed as a spin system, their geometry 
is very different from that of a traditional spin system on the lattice; our paper can be 
seen as a first step towards obtaining structural and mixing time results for triangulations 
analogous to those for classical spin systems. 

We mention finally that our mixing time result for the special case of 1-dimensional 
regions with A < 1 is related to work of Greenberg, Pascoe and Randall [7] on lattice 
paths. Those authors use a similar path coupling argument, but for a different probability 
distribution on lattice paths: in their model paths are biased according to the area under 
the path, while in ours the bias depends on the excursions of the path from a fixed line. 

2. The model 

2.1. Lattice triangulations. Let A^ nn denote the set of points in the mx n region of 
the integer lattice Z 2 , i.e., A^ n := {0,1,..., m} x {0,1,..., n}. A (full) triangulation 
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of A^ n is any maximal set of non-crossing edges (straight line segments), each of which 
connects two points of A^ and passes through no other point. We write n m>n for the 
set of triangulations of Aj^ . (We will drop the subscripts on A and f2 when their values 
are not important. Also, we will occasionally abuse notation by using A^ n to refer to the 
geometric region of M? that is the convex hull of the integer pointsO) In this section we 
spell out a few basic properties of Q m ,n that are either well known (see, e.g., [3]) or easily 
deduced. 

As for any point set in M 2 , the numbers of triangles and edges in every triangulation 
of A^ n are determined by the total number of points and the number of points on the 
boundary of the convex hull of the point set. Thus every triangulation in Q, m n has 
(m+ l)(n + l) points, 2mn triangles, and 3mn + m + n edges, 2(m + n) of which are on the 
boundary and the remainder in the interior. Moreover, since A^ n consists of the integer 

points inside a lattice polygon, every triangle is unimodular, i.e., has area |. 

In any triangulation, the midpoints of the edges are precisely the points of the half- 
integer lattice with the integer points removed, i.e., 

A mj n := {0, ^, 1, |, . . . , m — ^, m} x {0, ^, 1, |, . . . , n — ^, n} \ A^ n . 

Thus we may think of a triangulation as an assignment a = {cr x } X £\ m n of an edge to each 
point of A = A m ^ n . We call midpoints x with one integer and one half- integer coordinate 
"Type 1", and those with two half-integer coordinates "Type 2." Figure [3] shows some of 
the possible configurations for the edge at Type 1 and Type 2 points x. 





Figure 3. Possible configurations of the edge at x. The configurations 
differ slightly according to whether x is Type 1 (left-hand figure ) or Type 2 
(right-hand figure). For Type 1 points the minimum length configuration is 
vertical (or horizontal); for Type 2 it is a unit diagonal. 

Note that the minimum length configuration of the edge at a Type 1 point is horizontal 
or vertical, and at a Type 2 point it is a unit diagonal. The configurations of the boundary 
edges at points x = (0, j), (m,j) (respectively, x = (i,0), (i,n)) are forced to be horizontal 
(respectively, vertical) in all triangulations. We call triangulations in which all edges 
have minimal length ground state triangulations. There are exactly 2 mn ground state 
triangulations, which are equivalent up to flipping of unit diagonals. 

We will also consider a generalization to triangulations with the configurations of some 
interior edges fixed. Let A be an arbitrary subset of A, and rj = {^j^gA an arbitrary 
assignment of edges to the points of A that is consistent, i.e., none of the edges cross. 

^For notational convenience we shall adopt the slightly non-standard convention that denotes the 
point in R 2 with vertical coordinate i and horizontal coordinate j. Thus the region A^ n has vertical and 
horizontal dimensions m and n respectively. 
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We denote by A) the set of triangulations <r S O that agree with rj on A, i.e., all 
possible completions of the partial triangulation r/. Since triangulations are just maximal 
consistent sets of edges, it is clear that for any r\ at least one such completion always exists. 
We refer to the edges in r/ as constraints. See FigureHl 




FIGURE 4. Example of a triangulation of a 4x10 region. Above: constraint 
edges. Below: a triangulation consistent with the constraints (in bold). 

An important application of constraints is to triangulations of arbitrary lattice polygons 
(i.e., polygons — not necessarily convex — whose vertices are lattice points in Z 2 ). Let P 
be the set of integer points in a lattice polygon. Then triangulations of the point set P 
correspond to triangulations of an enclosing rectangle n with the polygon edges as 
constraints (and an arbitrary fixed triangulation outside the polygon). See Figure [5] for 
an example. 




Figure 5. Example of a triangulation of a lattice polygon embedded in a 
5x5 square. Left: the polygon edges (in bold), along with the edges of an 
arbitrary triangulation outside the polygon, are fixed as constraints. Right: 
a triangulation consistent with the constraints. 
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2.2. The flip graph. Any interior edge o x lies in two triangles. Prom area considerations, 
the third vertex of each of these triangles is an integer point on the line parallel to o x 
that passes through the closest integer point on either side of o x (see Figure 15(a)). The 
integer points on these lines occur periodically at intervals equal to the length of a x , and are 
positioned symmetrically on either side of o x , so that x is the midpoint of the line segment 
joining a point to its symmetric pair. The two triangles containing a x form a quadrilateral; 
this quadrilateral is convex iff the third vertices of the triangles are a symmetric pair. In 
this case a x is the diagonal of a parallelogram, and thus can be replaced by the other 
diagonal to yield another valid triangulation: o x is then said to be flippable. We observe 
that, for any a x that is not of minimum length, there is a a unique parallelogram in which 
o x is the longer diagonal; we call this the minimal •parallelogram of o x . Thus there is a 
unique flip that takes o x to a shorter edge. (There are many possible flips that make o x 
longer.) See Figure E(b,c). 




(a) (b) (c) 

Figure 6. Flips of an edge o x . (a) Possible locations of the third vertex 
of the triangles containing o x are the integer points on the dotted lines, 
(b) Unique flip of o x to a shorter edge (o x is the diagonal of its minimal 
parallelogram), (c) A flip of o x to a longer edge. 

We observe also that flipping an edge cannot change its slope from positive to negative 
(or vice versa), unless the edge is a unit diagonal. Thus any sequence of flips that changes 
the sign of the slope of an edge must pass through the minimum configuration of the edge 
(horizontal, vertical, or unit diagonal). Figure shows how this is achieved. 




Figure 7. Changing the sign of the slope of o x . (a) If x is Type 2, the 
change can be effected only when o x is a unit diagonal, (b) If x is Type 1, 
o x must first become horizontal (or vertical), and then be flipped after the 
orientation of its surrounding parallelogram has been changed. In both cases 
the new configuration is denoted o' x . 

The flipping operation induces a graph on vertex set 0, known as the "flip graph", in 
which two triangulations are adjacent iff they differ by one edge flip. The flip graph is 
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well known to be connected. To see this, note that a longest edge in any non-ground-state 
triangulation is always flippable to a shorter edge, because if a x is the longest edge in both 
of its triangles then the triangles must form the minimal parallelogram of a x . Hence, for 
any initial triangulation a, there exists a sequence of flips that reaches a ground state. 
Since flips are reversible, any triangulation is also reachable from a ground state. And 
finally, any ground state is reachable from any other by flipping unit diagonals. 

For triangulations with constraints we get a flip graph on the smaller set f2(r/, A). As we 
shall see in Section 13.21 in the presence of constraints there is still a well-defined ground 
state in which all configurations a x are of minimal length (subject to the constraints). 
Moreover, it can be shown (see Lemma 13.71 below) that any longest edge that is not in its 
ground state configuration is flippable to a shorter edge, and hence the flip graph remains 
connected in this case. In Proposition 13.81 we shall see that the flip graph actually enjoys 
a very strong structural property, which allows one to exactly compute the shortest path 
distance (or flip distance) between any pair of triangulations. 

2.3. Weighted triangulations and Glauber dynamics. We associate with a triangu- 
lation a a weight A'' 7 ', where A is a positive real parameter and |<r| denotes the total length 
of cr, i.e., \a\ = Y^ x \ a x\) here | • | denotes the i\ (Manhattan) metric. These weights induce 
a probability distribution \i on 0, (or, in the presence of constraints, on Q(rj, A)) via 

AM 

= y, (2-1) 

where the normalizing factor Z = ^ CT A' "' is the partition function. We refer to ()2.ip 
as the Gibbs distribution. The case A = 1 corresponds to the uniform distribution on 
triangulations; the cases A < 1 (respectively, A > 1) favor shorter (respectively, longer) 
triangulations. Note that when A < 1 the triangulations of maximum weight are precisely 
the ground state triangulations, i.e., those that minimize |cr|o 

The flip graph forms the basis of a local Markov chain (or Glauber dynamics) on f2 (or, 
more generally, on 0(r/, A)) as follows. In state a £ $7, pick a point x 6 A uniformly at 
random^; if the edge a x is flippable to edge a' x (producing a new triangulation a'), then 
flip it with probability 

/i(g') _ 
H(<r') + n(a) ~ AKI +AN' 

else do nothing. Since the flip graph is connected, this so-called "heat-bath" dynamics 
defines an ergodic Markov chain on Q (or on £l(n, A)) that is reversible with respect 
to \i. Hence the dynamics converges to the stationary distribution /j>. We will analyze 
convergence to stationarity via the standard notion of mixing time, defined by 

T mix = inf {t g N : max \\p\a, •) - fi\\ ^ 1/4} , (2.2) 

where ^((j, •) denotes the distribution after t steps when the initial state is cr, and /i|| = 
2 So-gc IK ") ~~ M ")! i s the usual total variation distance between two distributions /U, v. 



This explains our choice of the term "ground state"; in statistical physics, ground states are configura- 
tions of minimum energy, and thus maximum weight in the Gibbs distribution. The term "ground state" 
is less appropriate when A > 1, but our focus in this paper is mainly on the case A < 1. 

■^In the presence of constraints, we pick x u.a.r. from the midpoints of non- constraint edges. 
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3. Structural properties 

In this section we establish some basic geometric properties of lattice triangulations and 
the flip graph. Throughout we will work with arbitrary constraints, so that our results 
apply in particular to triangulations of any lattice polygon. 

3.1. The minimal parallelogram and excluded region. We begin with a useful prop- 
erty of the lattice Z 2 . 

Proposition 3.1. Let a,b be positive and coprime, and consider the infinite line through 
(0,0) of slope (a, b). Then the points of 7? that are closest to this line (and are not on the 
line) are at horizontal distance a -1 and vertical distance b^ 1 on either side of the line. 

Proof. Consider two successive integer points on the line (say, (0,0) and (a, &)), and the 
a horizontal lines at integer heights 0, 1, 2, . . . , a — 1. On each such horizontal line, con- 
sider the closest integer point to the line that lies strictly to the right of the line. (By 
symmetry, the points on the left of the line behave similarly.) It is readily verified that 
the horizontal distances of these points from the line are all distinct elements of the set 
{^:i = l,2,...,a}. Hence the closest of these points is at horizontal distance a -1 . This 
pattern repeats with period a throughout 1? . The same argument applies to the vertical 
distances, with period b. Finally, note that a point at horizontal distance a -1 is also at 
vertical distance b^ 1 . □ 

An immediate corollary is the following. 

Proposition 3.2. Let a, b be positive and coprime. The infinite family of parallel lines with 
slope (a, b) and horizontal separation a -1 (and thus vertical separation b" 1 ) partitions 1? , 
in the sense that every integer point lies on one of the lines. 

Given any triangulation a, consider an arbitrary edge o x that is not horizontal or 
vertical. Let its slope be (a, b) with a, b coprime, and denote its upper and lower endpoints 
by Pi , Pi respectively. By Proposition 13.11 and its proof we may identify the unique closest 
integer point to the right of a x that lies horizontally and vertically between p\ and p2- 
Call this point p^. By symmetry there is a corresponding closest point p^ on the other 
side of a x , so that x is the midpoint of the line segment p^Pa- The points Px,P3,P2,P4 form 
the vertices of the minimal parallelogram of a x mentioned earlier: o x can be flipped to the 
shorter edge p^P^ iff all edges of the minimal parallelogram are present in the triangulation. 

The following statement identifies a region centered on the minimal parallelogram that 
is free of integer points. (See Figure 0) 

Proposition 3.3. Let a x be any non-horizontal and non-vertical triangulation edge with 
endpoints pi,P2 and minimal parallelogram p\p^P2PA- Then the shaded region in Figured 
(the "excluded region") does not contain any integer points in its interior. 

Note: The excluded region is the union of two strips with parallel sides that extend to the 
boundaries of A . The intersection of the two strips is the minimal parallelogram of a x . 
The complement of the excluded region consists of four disjoint closed cones, which we 
shall call R\, R2, R3, R4, as shown in Figure M (so the apex of Ri is Pi). 

Proof. It suffices to show that each of the two pairs of parallel lines that bound the strips 
are adjacent lines in a family as in Proposition 13.21 Consider first the pair of lines through 
P1P4 and P2P3- We need to verify that p\ is a closest point to the line through P2P3- 
W.l.o.g., assume that P2 is the point (0,0) and p\ is the point (a, b), with a, b positive 
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Figure 8. The minimal parallelogram and excluded region of an edge o~ x . 

and coprime. Suppose P3 is the point (a', b'), and recall that the horizontal distance of p^ 
from a x is a -1 . By similarity of triangles, the horizontal distance of p\ from the line 
through (distance p\u in Figure [8]) is x - = \. Hence, by Proposition 13.11 p\ is 
a closest point to this line, and by Proposition 13.21 this parallel strip contains no integer 
points. 

An analogous argument using vertical separation shows that the distance P2V is ^ and 
hence the other parallel strip, bounded by lines through P2PA and pip?,, also contains no 
integer points. This completes the proof. □ 

3.2. The ground state lemma. We turn now to the structure of the ground state trian- 
gulations, i.e., those of minimum total length. These are the triangulations of maximum 
weight in the probability distribution (|2.ip when A < 1, and they play a central role in 
our analysis of spatial mixing in Section [H In the absence of constraints, the ground 
state triangulations are very simple: every edge is either horizontal or vertical or a unit 
diagonal, so in particular the ground state is unique up to flipping of the unit diagonals. 
The presence of constraint edges, however, may change the ground state considerably. In 
this subsection we prove that the ground state remains (essentially) unique, and can be 
described easily (see Figure [9] for an example). 

Lemma 3.4. [Ground State Lemma] Given any set of constraint edges, the ground 
state triangulation is unique (up to possible flipping of unit diagonals), and can be con- 
structed by placing each edge in its minimal length configuration consistent with the con- 
straints, independent of the other edges. 

Lemma 13.41 will follow from a short sequence of observations. Given a set of con- 
straints (77, A), we call a midpoint x constrained if the minimal configuration of the edge 
at x consistent with the constraints is not horizontal or vertical or unit diagonal; or if the 
minimal configuration is a unit diagonal and its other unit diagonal configuration is not 
consistent. 

Let a x be a consistent configuration for the edge at x with endpoints pi,f>2 as in Figure [HJ 
We say that a x is spanned by a constraint edge C if C crosses the line segment P3P4 in 
Figure [8] (i.e., the shorter diagonal of the minimal parallelogram of a x ). Note that this 
necessarily implies that the endpoints of C lie in the regions R\ and i?2- 
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Figure 9. Example of a ground state triangulation in the presence of the 
constraints given in Figure^ 

Proposition 3.5. Given a set of constraints, a consistent configuration o~ x for the edge 
at a constrained x is minimal if and only if a x is spanned by a constraint edge. Moreover, 
the minimal configuration a x = a x is unique. 

Proof. Assume first that a x is minimal and x is constrained. Then in particular the edge 
P3P4 must not be a consistent configuration for the edge at x (since either it would be a 
shorter configuration, or in the case of a unit diagonal it is explicitly prohibited). Hence 
a constraint edge must pass through P3P4. 

For the other direction, suppose a x is spanned by a constraint edge. We claim that 
any consistent configuration for the edge at x must have its endpoints in the cones R\,R-2 
defined by a x (see Figure E]). This immediately implies that a x is the unique minimal 
configuration for this edge. To see the above claim, let C be the spanning constraint edge. 
Since any configuration a x for the edge at x must be obtained from a x by rotating about 
its midpoint x, and there are no integer points in the excluded region, the endpoints of a' x 
must lie either in R\ , R2 or in Rj,,R^. But an edge between R3 and R4 would necessarily 
cross C, so we are left only with the possibility R\,R2- □ 

Proposition 3.6. Given a set of constraints, let o x ,a y be minimal configurations of any 
two distinct edges. Then cr x ,a y do not cross. 

Proof. By Proposition 13.51 both a x and a y are either spanned by a constraint edge or the 
corresponding midpoints are unconstrained. If both x, y are unconstrained then clearly 
they do not cross. So assume w.l.o.g. that a x is spanned by some constraint edge C. 
Suppose for contradiction that o~y crosses o~ x . Then o~y must have endpoints either in 
R\, i?2 or in .R3, i?4, where the Ri are the regions defined by o x as in Figure As in the 
proof of Proposition 13.51 an edge between .R3 and R4 would necessarily cross C, so a y 
must have endpoints in -Ri , i?2 - But clearly this implies that \a y \ > \a x \. 

Now if y is constrained, then switching the roles of x and y in the above argument yields 
\&x\ > l^i/l) a contradiction. And if y is unconstrained then a y is horizontal, vertical or 
a unit diagonal, in which case it cannot be strictly longer than the constrained edge a x , 
again a contradiction. □ 

The Ground State Lemma (Lemma l3.4p now follows immediately from Propositions 13.51 
and 13.61 Given any set of constraints (77, A), we may thus speak of a minimal length 
triangulation a = {cr x } X £\ which is unique except for possible flipping of unit diagonals. 
We call a a ground state triangulation. 

3.3. Flip distance. We start with a proof of the claim in Section \2. 21 that the flip graph 
remains connected in the presence of arbitrary constraints. (Recall that connectedness is 
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easily verified in the absence of constraints.) We shall also see that the diameter is always 
small, regardless of the constraints. 

Lemma 3.7. For any set of constraints (rj, A), the flip graph on triangulations in Cl(r], A) 
is connected. Moreover, its diameter is 0{mn{m + n)). 

Proof. By the Ground State Lemma (Lemma I3.4p , there is a well-defined family of ground 
state triangulations for any such (rj, A), which differ only up to flipping of unit diagonals. 
As in the unconstrained case, it suffices to show that a ground state is reachable from 
any triangulation in fi(ry, A), since by reversibility this implies that any triangulation is 
reachable from a ground state, and ground states are obviously reachable from each other. 
Since in any ground state triangulation all edges have minimal possible length, it suffices in 
turn to prove that, in any non-ground-state triangulation, there is an edge that is flippable 
to a shorter edge. This we now do. 

Let a G £2(77, A) be an arbitrary non-ground-state triangulation, and let a x be some edge 
of a that is not in ground state. If a x is flippable to a shorter edge we are done; if not, 
then o x cannot be the longest edge in both of its triangles, so we can find a longer edge, 
o~ y , that shares an endpoint with a x . The crucial observation (see the next paragraph) 
is that o~ y must also not be in ground state. Thus we can iterate this process, finding a 
connected sequence of edges of increasing length that are not in ground state, until we 
find one that is longest in both its triangles. This edge must be flippable to a shorter edge 
and we are done. 

To conclude, we verify the above claim that a y must not be in ground state. This 
follows from the fact that, in any triangle, if the longest edge is in ground state then so 
are the other two edges. To see this, simply note that the longest edge in a triangle spans 
the other two edges; this can be seen because the triangle contains no lattice points, so 
the longest edge must pass between each of the other edges and its closest lattice point on 
the same side of the edge (see Figure [8]) . 

Finally, the assertion concerning the diameter follows immediately from the fact that 
the maximum length of a triangulation is 0(mn(m + n)), which is clearly an upper bound 
on the number of edge-shortening flips needed to reach a ground state triangulation. □ 

Next, we shall compute the flip distance between two arbitrary triangulations in the 
flip graph, i.e., the minimal number of flips required to obtain one triangulation from the 
other. As usual, fix a set of constraints (77, A). Given a midpoint x, let £l x denote the 
set of all possible values of the edge at x that are consistent with the constraints. This 
set contains the ground state o x (or, in the case where x is unconstrained and of Type 2, 
the pair of ground states corresponding to two opposite unit diagonals). Say that two 
elements a x , a' x G £l x are neighbors, written a x ~ o~' x , if a x is flippable to a' x within a valid 
triangulation a G £2(77, A). Recall that for any a x G £l x there is at most one o~' x ~ a x 
such that \o~ x \ ^ \o~' x \. If a x is not in ground state, then an edge a' x G Sl x with o' x ~ o x 
and \o x \ > \a' x \ necessarily exists. This follows from the process described in the proof of 
Lemma 13.71 since o x must eventually be flippable to a shorter edge. These observations 
prove that the resulting graph with vertex set Vl x is a tree rooted at the ground state o x (or, 
when the ground state is not unique, at a pair of neighboring ground states corresponding 
to unit diagonals). Given o~ x ,t x G £l x , let k(o~ x ,t x ) be the distance between o~ x ,t x in the 
tree, i.e., the minimal number of flips of the edge at x required to change its configuration 
from o x to t x (regardless of the disposition of the other edges). Notice that any path 
between a and r in the flip graph must contain all of the above flips for every x, and thus 
have length at least Yl x k (o~x,t x ). We call these the indispensable flips for the pair a, r. 
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The next result expresses a fundamental structural property of the flip graph: there is 
a path between any two triangulations that consists only of indispensable flips. Note that 
this path always exists, independent of the constraints. Although we shall not require this 
result in our subsequent analysis, we include it here for completeness. 

Proposition 3.8. The flip distance between any two triangulations o~,t £ Q(r/, A) is equal 

Proof. Given any pair of triangulations o~,t G A) with a ^ r, we show how to 
perform an indispensable flip in either a or r to produce a new pair o~',t'. Since the set of 
indispensable flips for a', r' differs from that for a, r only in the single flip just performed, 
iteration of this procedure until a = r produces a path of indispensable flips between a 
and t. 

To justify the above claim, let A' = {x E A : a x = t x }; note that A' contains A, and 
(A, rj) = (A, a). We will view (A',ct) as an expanded set of constraints. Let denote a 
point x G A \ A' such that 

\a X)t \ = max-fjo-a;! : xGA\A'}. 

Similarly, let denote a maximal edge in r outside A'. Assume w.l.o.g. that \cr x J > \r yt \ 
(else the same argument follows with the roles of a Xit and r yt interchanged). We claim 
that an indispensable flip can be performed on o~ Xt . 

Note first that \a Xt \ > \t v J > \t x J, so since ground states are minimal and o Xf / r Xtt , 
the only way a Xf can be in ground state is if |<r x J = \t x J and a Xt ,T Xt are opposite unit 
diagonals; in this case an indispensable flip of a Xt must be possible (since it could be 
prevented only by a longer non-constraint edge, and a Xf is assumed to be maximal) and 
takes us directly to T Xt . Otherwise, o~ Xt is not in ground state, so from the proof of 
Lemma 13.71 since a Xf is a longest such edge it is flippable to a shorter edge. To see that 
this flip is indispensable, note that a Xt / T Xt and \a Xt \ > \r Xit \, so the path from a Xt to 
r Xf in the tree £l x must include the shortening flip at o~ Xt . This completes the proof. □ 

3.4. The influence region. The following notion of "influence region" of an edge will 
play a key role in our analysis of decay of correlations under the Gibbs distribution for 
small A. As usual we fix an arbitrary set of constraints (r/, A). Let ex be a ground state 
triangulation compatible with the constraints as in Lemma l3.4i Suppose we change the 
configuration of the edge at some point x to a non-ground-state value a x ^ a x (still 
consistent with the constraints). We would like to identify the minimal region "affected 
by" the new edge a x . We call this region the influence region of a x , and denote it T(a x ). 
In this section we give two definitions of T(a x ) — one implicit, the other constructive — and 
prove their equivalence. 

For the first definition, we define T(a x ) as the smallest region S C M? made up of 
ground state triangles such that a x is contained in S. (When a x = a x , we define T(<r x ) to 
be empty. Note that the ground state triangles in S may differ from those in a up to flips 
of unit diagonals.) Observe that, since the boundary of T(a x ) must consist of ground state 
edges that do not cross a x , the presence of a x has no effect on the ground state outside 
this region. 

Now we describe an equivalent construction of the set T(a x ) (refer to Figure [T0|) . In the 
setting in the first paragraph above, let a(a x ) denote the new ground state triangulation 
when the edge a x is added as an additional constraint. Consider the region T*(a x ) C M 2 
defined by successive addition of triangles as follows. Let yi,y2 and Z\,Z2 denote the 
midpoints of the edges of the minimal parallelogram of o~ x , so that this parallelogram is 
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Figure 10. The influence region T(cr x ) when o~ x is an edge of slope (—3, 5) 
(with midpoint in the center of the box). Left: in the case of no constraints. 
Right: in the presence of a constraint edge of slope (—4,7). Above: the 
triangles inside correspond to the second definition of the influence region. 
Below: the triangles inside correspond to the first definition of the influence 
region. In the background, a ground state triangulation (dotted lines). 



made up of the two triangles T\ = {a x ,a yi (a x ),a y2 (a x )) and T 2 = {(7 x ,a Zl (a x ),a Z2 (a x )). 
Note that these triangles are consistent with the constraints: indeed, if one of the con- 
straint edges C were to cross an edge of one of these triangles, then C would span a x ; 
and by Proposition 13.51 this would imply that a x is minimal, contradicting the assump- 
tion that a x ^ a x . Thus, we can add the triangles T\ and T 2 . Next, we add a further 
triangle for each edge a Vt (a x ), i = 1,2, such that a Vi (a x ) ^ a Vi and for each edge a Zi (a x ), 
i = 1,2, such that a z% (a x ) ^ a Zi . For edges that are already in their original ground 
state we add no further triangle. More precisely, if a yi (a x ) ^ a yi , then define mid- 
points yii, 1/12 and the triangle T%\ = (a yi (cr x ), cr yn (cr x ), cr yi2 (cr x )), which is one half of 
the minimal parallelogram of a yi (a x ). The same reasoning as above shows that this tri- 
angle is consistent with the constraints as long as cr yi {a x ) ^ a yx . Similarly, define 2/21,2/22 
and T12 = (^ y2 (a x ),a y21 (a x ),a y22 (a x )), which is one half of the minimal parallelogram of 
a y2 (o~ x ). The same construction can be applied to the other side of a x to define the points 
0U, z\2, Z21, Z22 and the associated triangles T21, T22. This branching procedure is iterated 
until one of the midpoints, say w, is such that cr w (a x ) = a w . In this case that branch is 
stopped and we continue with the other available branches, if any. T*(a x ) is defined as 
the union of all triangles added in this way. 

Notice that all triangles added contain at least one edge that is not in ground state. 
However, the boundary of Y*(o~ x ) is made up only of ground state edges. Therefore, 
T*(a x ) is also a union of ground state triangles. (See Figure [TUl) 

Lemma 3.9. For any point x and any value of the edge a x , Y*(o~ x ) = T(o~ x ). 

Proof. We can assume o x 7^ a x , since otherwise T*(a x ) = T(a x ) = 0. Since T*(a x ) is a 
union of ground state triangles that contains a x , we have T(a x ) C Y*(o~ x ). Next, suppose 
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for contradiction that T(a x ) is strictly contained in T*(o~ x ). Then there must be a point z 
such that a z is a boundary edge of T(a x ) but z is an interior point of T*(a x ). This implies 
that z can achieve its ground state value a z even in the presence of a x , but the fact that z 
is an interior point of T*(a x ) implies that a z {a x ) 7^ a z (otherwise the branching procedure 
defining T*(a x ) would have stopped at z and z would not be an interior point). Thus, 
such a z cannot exist. □ 



4. Spatial mixing for small A 

In this section we prove some fundamental properties of the Gibbs distribution fi for 
sufficiently small values A < 1. Since long edges are penalized in this regime, one might 
expect that the probability that an edge is longer than its ground state value should 
decrease exponentially in this excess length. We prove this property in Corollary 14.21 
below, as a consequence of a rather more general exponential tail property (Lemma 14. ip . 
We then go on to establish exponential decay of correlations (or spatial mixing) in the 
same regime (see Theorem 14. 4p . as advertised in Theorem 11.21 in the Introduction. It is 
worth emphasizing that all these properties are shown to hold uniformly in the constraints 
(77, A), and thus in particular for triangulations of arbitrary lattice polygons. 

We begin with some additional terminology. Fix a set of constraint edges. For any 
triangulation a consistent with the constraints, we call a triangle in a a G-triangle if all of 
its edges are in ground state (with respect to the constraints), and a B-triangle otherwise. 
Two triangles are connected if they share an edge, and ^-connected if they share a vertex. 
Note that every maximal connected region of B-triangles is bounded by ground state 
edges (see Figure fTT]) . For any midpoint x, we denote by S x (a) the maximal connected 



FIGURE 11. The triangulation with constraints of Figure^ showing its 
four connected regions of B-triangles (shaded). All unshaded triangles are 
G-triangles. Compare also with Figure^ 

component of B-triangles containing the edge at x (so S x (o~) = iff a x is contained in two 
G-triangles). Notice that S x (o~) need not be simply connected (i.e., there may be holes 
made up of G-triangles), and the boundary of S x (a) consists of ground state edges. We 
call S = S(o~) the union of all B-triangles in a, i.e., S(a) = \J X S x (a). 

For any region SCA° that is the union of ground state triangles, let A(S') denote the 
set of interior midpoints of S, namely the collection of z G A n S that are not on the 
boundary of S. Notice that any triangulation of the region S is obtained by specifying 
the edges a z , z £ A(S). We define 

*(S,o):= 1^1- 

zeA(S) 
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In particular, $(S x (o~),o~) is just the total length of the triangulation inside the region 
S x (a). Clearly, for any a, we have $>(S x (a),a) ^ $>(S x (a), a), where a denotes a ground 
state triangulation. The inequality is strict whenever S x (cr) 7^ 0. It will be useful to 
consider the quantity 

</> x (a) :=${S x (a),a)-${S x ((T),a), (4.3) 
The next lemma shows that, for sufficiently small A, 4> x has exponentially decaying tail 
probability. 

Lemma 4.1. There exists Ao £ (0, 1) such that, for all sets of constraints, all points 
x e A, all A > and all k G N, 

fi(M°) = k)^(\/\o) k - (4.4) 

Proof. We have 

n((j) x (o-) = k) = } ] jj,((f> x (o-) = k;S x (a) = S). 
s 

Let Z denote the partition function, i.e., Z = X^o-era ^' <T ' • ^ or a fi xe d region S that is 
the union of ground state triangles, write r ~ S if r is a triangulation of S, and let 
Zs = X^r~5 ^' T ' be- the restricted partition function within S. As a convention, for r ~ S, 
we do not include the boundary edges in r, so that 

\t\=<S>(S,t)= \tz\. 

zeA(S) 

Notice that for any region S as above we have the bound 

Z > \\ ds \ZsZ S c, 

where \dS\ denotes the length of the boundary of S. For a triangulation a ~ S we also 
write <t~*S' if every triangle in a ~ S has at least one edge not in ground state. Then, 



z 

z 

Y^a^S ^ <I>( '' S '' CT ^l'I»(5,o-) = *(5,5-)+fe; cr~*S 



Z s 

Next, we can use the trivial bound Z$ ^ \®( s >°j } so that 

S <7~S 

Since a ~* S, we have $>(S,a) — <fr(S,a) ^ ^\S\, where |5| stands for the number of 
triangles in S. (Recall that, since all triangles are unimodular, this number does not 
depend on the particular triangulation a ~ S). Then, for any Ao > we have 

V(<p x (a) = *) < (A/A ) fe £ £ Ajf l/2 . (4.5) 

S (T~S 

From Anclin's bound [1] mentioned in Section Tl .31 there are at most 2 3 ' ,s ' //2 triangulations 
of a region S, since the number of interior edges is at most 3|S'|/2. Moreover, the number 
of connected regions S containing the point x as above and with a given size \S\ is bounded 
by9l 5 l. To see this, observe that every such S must be a collection of \S\ connected ground 
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state triangles, so it suffices to count the number of connected sets made up of i triangles 
in a valid triangulation, which is at most 3 2 ^ since each triangle is connected to at most 3 
other triangles and a general connected set can be explored by a path that visits the same 
point at most twice. In conclusion, the claim follows by summing in (|4.5h : 

oo 

KM*) = k) < (X/Xo) k £ 2^X[ S ^ < (A/A ) fe CX /2 , (4-6) 

s e=2 

where Co = 9\/8 and we use the fact that the minimal non trivial S has size at least 2. It 

1/2 

suffices to take Ao > such that (e.g.) A Co ^ 1/2 to ensure that the last summation is 
less than 1. □ 

We are now able to prove the exponential tail bound advertised at the start of this 
section. Notice that if a G f2 is such that |<r z | = \a x \ + k, for some point x G A, then 
necessarily (j) x (a) ^ k. Hence from Lemma 14.11 we immediately have: 

Corollary 4.2. Let Xq be as in Lemma \4.1\ For every point x £ A and every k £ N, 

KWxl = \° x \ + k) ^ (A/A ) fc . (4.7) 

Recall that S = S(a) is the random subset of the plane defined as the union of all 
B-triangles in a. The set S naturally describes the deviation of a from the ground state a. 
Our next result, based on a slight variation of the argument of Lemma [4. 11 establishes that 
for small A the probability that S contains any particular region V decays exponentially 
in the size of V. Since the boundary of S consists of ground state edges, S can always be 
described as a collection of ground state triangles. Thus we naturally also specify V as a 
collection of ground state triangles. We write \V\ for the number of triangles in V. 

Lemma 4.3. Let Ao be as in Lemma \4-l\ For every V C A that is the union of ground 
state triangles, we have 

KV C S) < (4A/A ) |y|/2 . (4.8) 

Proof. Suppose first that V is connected. Let T be a ground state triangle in V and let 
So denote the connected component of S containing T. Then 

s-.sdv 

where the sum is over all connected sets S of ground state triangles such that S 5 V. 
Reasoning as in the proof of Lemma 14.11 we have that 

^(S = 5) < (2 3 A) |5|/2 . (4.9) 

Summing over all S as above and estimating by 9^ the number of connected S B T with 
\S\ = £, with the notation of (|4.6p we have 

Y C^/^(X/X r^. 
t>\v\ 

We turn to the case where V has several connected components, say V\, . . . ,V n . Let 
Ti,...,T n denote fixed ground state triangles such that Tj G Vi, and let Si, i = 1, . . . , n 
denote the random sets defined as the connected components of S such that Si 3 Tj. 
Notice that the Si need not be distinct, since a single connected component of S may 
contain more than one of the triangles Tj. However, the event V C S implies that there 
exist i%,...,i m , 1 ^ m ^ n, and distinct connected components S^, . . . , Si m such that 
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Si k 3 Ti k and Si k = Si k , for all k = l,...,m. For any 77 G such that Si k = Si k , 
k = 1, . . . , m, reasoning as in the proof of Lemma 14,11 and (|4.9p we have 

fi(S ik =S ik , Vfc = l,...,m) < (2 3 A)^™=il^l/ 2 . 

1/2 

Therefore, using CoAp 7 < 1/2, n ^ |V|, and a union bound we arrive at 



m=l v 7 ^ > 1 £ m > 1 L J 



(CoA i/2)Er =1 4 



< E (™V/ A °) |V1/2 < 2n (VAo) |y|/2 < (4A/A )l y l/ 2 . 

m=l 

□ 

Finally, we use Lemma 14.31 to derive the spatial mixing estimate stated in Theorem 11,21 
in the Introduction. This bound will be expressed in terms of a natural distance between 
the edge a x and any subset of midpoints A C A, which we now define. 

Let a x be an edge consistent with the constraints (rj, A). Also, let a(a x ) denote the new 
ground state triangulation when edge a x is added to the constraints. Let fi ax denote the 
Gibbs distribution on triangulations with this extra constraint, i.e., (i ax is the probability 
H conditioned on the event that the edge at x equals a x . Intuitively, the distributions pf x 
and [i should have very similar marginals at midpoints z that are far from the influence 
region T(a x ), as defined in Section 13.41 To quantify this we introduce the following 
distance. Consider the graph Q{a x ) whose vertices are the triangles in a(a x ), where two 
triangles are neighbors iff they share an edge. Let d(A,o~ x ) denote the graph distance in 
G(a x ) between the influence region T(a x ) and the set of triangles spanned by edges z 6 A 
(i.e., the union of all triangles in a(o x ) containing the edges a z (a x ), z £ A), where the 
distance between two sets of vertices of a graph is interpreted as the minimum distance 
between a pair of vertices, one from each set. In words, d(A,a x ) is the minimal number 
of triangles in cr(a x ) needed to connect an edge in A to the influence region T(a x ). Below 
we use I \fi r7x — hWa for the variation distance between the distributions fi Ux and fi at A, 
that is the variation distance between the marginals on edges z £ A. We write \A\ for the 
cardinality of A. 

Theorem 4.4. There exists Ai € (0, 1) such that, for all sets of constraints, all x G A 
and AC A and all consistent values of the edge a x , we have 

||/i ff --/i[U<|^|(A/A 1 )^-). (4.10) 

Proof. Observe that apart from the set T(o~ x ) the ground states for \x ax and fi coincide. 
When we talk about triangles below we always refer to the triangles in the common ground 
state triangulation (outside the region T(a x )). Moreover, we identify A with the union of 
all triangles spanned by edges from A. Let P denote the independent coupling of \i ax , \i 
and let a, a' denote the associated random triangulations. Consider the event £ that there 
exists a *-connected chain C of ground state triangles that separates T(a x ) from A such 
that both a and a' contain C. From the Markov property of the Gibbs distributions \i ax , \i 
it follows that they can be coupled at A in such a way that 

IK* - fi\\ A ^ F(£ c ). (4.11) 

Next, observe that the complementary event £ c implies that there exists a connected 
chain T> of triangles connecting T(a x ) and A such that, for every triangle T G D, either 
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T £ S(a) or T G S(a'), where S(cr) as usual denotes the set of B-triangles in a triangulation 
a. Thus, given the set T>, either <S(cr) or S(o~') contains at least half of the triangles in D. 
By definition \D\ ^ d(A,a x ). In particular, from the union bound it follows that 

IK*-HU<E E \m°'(Yqs) + h(vqs)] 

V VCD: \V\ > \V\/2 

where the first sum ranges over all connected subsets T> touching A and such that \T>\ d(A, a x ). 
From Lemma S3] it follows that for any set V, n ax (V C <S) < (4A/A ) |y|/2 and //(V C 
5) < (4A/A ) |y|/2 , and therefore 

E \M (T '(VCS) + t ji(VQS)] ^ 2l^l(4A/Ao) |I?l/4 . 

VCV:\V\ S> \V\/2 

Using the fact that the number of connected T> touching a given triangle T C A with 
|P| = k is at most 9 fc , and bounding by \A\ the number of T C A, we get 

IK- -MlU < 1^1 E ^(4A/Ao) fc/4 ^ |A|(A/Ai) d ^)/ 4 , 

for suitable constants Ci,Ai > 0. □ 

5. Upper bounds on the mixing time for A < 1 

We consider the Glauber dynamics on triangulations of the region A := n , as defined 
in Section \2. 31 We allow an arbitrary set of constraints (rj, A), so that our results apply in 
particular to triangulations of any lattice polygon. Our main result in this section is the 
following polynomial upper bound on the mixing time for all sufficiently small A, which 
was stated as Theorem II .31 in the Introduction. (Recall the definition of mixing time from 
equation (12.2p .) 

Theorem 5.1. There exist constants C > and Ao € (0, 1) such that, for all A ^ Ao, all 

m,n G N, and all constraints (rj, A), 

T m ix ^ Cmn(m + n). (5-12) 

We first establish a key coupling estimate and then turn to the proof of Theorem 15.11 

5.1. Path Coupling. Let a and r be two triangulations so that a differs from r only on 
the edge with midpoint x. Let a > 1 be a fixed positive number. If a x and t x are opposite 
unit diagonals, define A(<r, r) = a 2 — 1; otherwise, let 

A(o-,t) = \a M - a^\. (5.13) 

Note that A is defined only over pairs of triangulations that are adjacent in the flip graph. 
We extend A to all pairs of triangulations by defining A(<r, r) to be the shortest path 
distance between a and r in the flip graph, with edge lengths given by the above A values. 

Lemma 5.2. There exists ao > 1 such that the following holds for all a ^ ao, and 
Aa ^ 1. Let a,r be two triangulations that differ in exactly one edge and let o~',t' denote 
the random triangulation obtained after one step of the Markov chain. Then, there exists 
a coupling of a' and r' so that their expected distance satisfies 

E[A(^rOKA(a,r)(l-^), (5.14) 

where A is the set of midpoints of non- constraint edges. 
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Proof. Let x be the midpoint so that o x and r x differ. Assume without loss of generality 
that |er x | ^ \t x \. Since a differs from r by one flip at x, the edges a x and t x must be 
the two diagonals of the minimal parallelogram of a x . The coupling is such that the 
random midpoint Y £ A chosen is the same for both triangulations. Thus, if Y = x then 
a' x = t' x and A(<r / ,r / ) = with probability one. If instead F/x, which happens with 
probability (1 — 1/|A|), then A(a' , t') = A(er, r) unless Y is one of the edges of the minimal 
parallelogram of a x (since only these edges share a triangle with the differing edge at x); 
so it remains only to bound the expected change in distance induced by flips at the edges 
of the minimal parallelogram. 

Assume first that \a x \ > \t x \. Note that in this case the length of a x can be written as 
fa + fa , where fa and fa are the lengths of the other edges of the triangles containing a x 
in a. Assume w.l.o.g. that fa ^ fa; then we have \t x \ = fa — fa- Let y\ , 1/2 , z\ , z<i G A be 
the midpoints of the edges of the minimal parallelogram, so that \a yi \ = \a y2 \ = fa ^ fa = 
Wzi \ = \ <J z 2 \- We start with the case Y = z\. Note that if a zi (respectively, T Z1 ) is flippable, 
then t Zi (respectively, a Zl ) is not flippable. Since fa ^ fa, if cr Zl is flippable it must flip 
to an edge a' zi of length \a x \ + \a yi \ = 2 fa + fa- Then, for a > 1, A < 1 with aX ^ 1, the 
expected increase in distance in this case is bounded by 

( a 2/ 1+ / 2 _ a / 2 ^ A2/l+/2 _ a^(l-a-^)(Aa)^ < ^ 
V ° J (\2A+h+\h) (A 2 /i +1) ^ • 

If instead t zi is flippable, then it flips to an edge of length \t x \ + \r yi \ = 2 fa — fa, resulting 
in an expected increase in distance of at most 

Now we turn to the case Y = y\ . Note that a yi is either unflippable or flippable to a longer 
edge, while r yi is either unflippable or flippable to a shorter edge. If o~ yi is flippable, it 
flips to an edge of length \a x \ + \o~ zl \ = fi + 2 fa, and the expected increase in distance is 
at most 

\/l+2/ 2 

If T yi is flippable, it flips to an edge of length \t x \ + \t Zi \ = \fa — 2fa\, which causes an 
expected increase in distance of at most 

^ a l/i-2/a| _ a h^ ^ a l/i- 2 /2l ^ a h . 

Summing up all these contributions, we obtain 

„.,,, hl A , J 1 \ 2ah 2a h 



This implies the desired bound in the lemma, since if a is large enough (a > 5 suffices) 
we have 

2a h + 2a h ^ - {a h+h - a h ~ h ) = - A(cj, r), 

for any fa > fa ^ 1. 

It remains to consider the case |<r x | = \t x \; note that here A(a,r) = a 2 — I. This 
case can only happen if o x ,t x are opposite unit diagonals, so \o x \ = \t x \ = fa + fa with 
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fl = fi = 1. Therefore, the total contribution from the edges of the minimal parallelogram 
(now a square) is bounded by 

A 3 

4 (a 3 -a) ^ 3 - A ^ < 4a -1 (a 2 - l)(Aa) 2 < 4a" 1 A(<r, r). 
This gives the following bound on the expected distance after one flip: 
E[A(a'y)]„A(a,r)(l-l J ) + Aa - lA( ^ ) 



\MJ |A| ■ 

The desired bound again follows by taking a ^ 8. □ 

Note: In the above proof it suffices to take a = 8 and A ^ a -1 = 1/8. Hence Theorem 15. II 
holds with Ao = 1/8. 

5.2. Proof of Theorem 15.11 Once the estimate in Lemma [5. 21 is available the argument 
is rather standard; see, e.g., [11} Theorem 14.6]. Indeed, from the triangle inequality and 
the definition of the metric A(-,-), it follows that the estimate (I5.14D can be extended 
to any pair of triangulations a, r under a suitable coupling. This fact, together with the 
Markov property yields the bound 

E[A(aW,rW)] ^ A(a,r) (l-J-) ^ A(er, r) e~ fc / 2 ' A ', (5.15) 



2|A| 

where E denotes expectation over the coupling of the random triangulations 
obtained after k steps of the Markov chain started at a, r respectively. The diameter bound 
in Lemma 13.71 and the fact that the distance between any two adjacent triangulations is 
at most a m+n , imply that A(<7, r) < e c °( m+n * > for some constant Co > 0. Using also the 
fact that A(<7, t) ^ c if a 7^ r, for some constant c = c(a) > 0, an application of Markov's 
inequality yields 

\\p k {a r ) -p fc (T,-)|| ^ c -l e O)(m+n)-fc/2]A] ; 

uniformly in the initial conditions a, t. Taking, e.g., k = 4|A|Co(m + n) concludes the 
proof. □ 

Note: Observe from the above proof that Theorem 15.11 actually holds with the bound 
Imix ^ C|A|(m + n), where |A| is the number of non- constraint edges. This bound may 
be much better in cases where there are many constraints. 

5.3. The 1 x n case for all A < 1. Here we prove that, in the special case of the 1- 
dimensional region A® n , the bound of Theorem 15.11 holds for all A < 1. 

Theorem 5.3. Set m = 1. For any A < 1, there exists a constant C > such that for all 
n E N, and for all constraints (77, A), 

T mix ^ Cn 2 . (5.16) 

To prove the theorem, we show that the path coupling argument of Section 15.11 works 
for all A < 1 in this case. For ID triangulations we may define the length of an edge as 
its horizontal length. Fix a > 1. Let a, r be two triangulations that differ at exactly one 
edge x with \a x \ ^ \t x \- Then either the edges at x have equal length, in which case they 
are necessarily both unit diagonals, or they have lengths I — 1 (in r) and I + 1 (in a), for 
some integer £ ^ 1. In the first case we set A(<r, r) = a 2 — 1, in the second case we set 
A(<t, t) = c/ +1 (l — a~ 2 ). As before, we extend A to the shortest-path metric on all pairs 
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in f2. As detailed in Section \5.2\ the estimate (|5.16|) will follow from the lemma below and 
the fact that the maximal distance between two arbitrary configurations is e°^ n \ 

Lemma 5.4. For any A < 1 there exist constants a > 1, 5 > such that the following 
holds. Let a, r be two triangulations that differ at exactly one edge, and let o~',t' denote 
the random triangulation obtained after one step. There exists a coupling of a',r' such 
that 

E [A(a', /)] < (1 - S/n)A{a, r) (5.17) 

Proof. Say that the initial discrepancy is at midpoint x. In the coupling, we pick the same 
random midpoint X to be updated, so that E[d(a' , r')] is given by 

E[A(cr', /)] = (1 - l/n)A(a, r) + - U(a, r), 

n 

where U(o~, r) is the term coming from possible new discrepancies at x ± 1/2. To compute 
the latter, we start with the case where o~ x ^ t x are both unit diagonals, so that A(<r, r) = 
a 2 — 1. In this case, the edges at x ± 1/2 are necessarily vertical in both er, r and we can 
create a discrepancy at either x ± 1/2 by letting the edge there increase to length 2. Only 
one of the two triangulations can flip at either x ± 1/2. Thus, taking expectations, in this 
case we have 

Next, consider the remaining case where \a x \ =1+1 and \t x \ = £ — 1. Necessarily the 
edges at x ± 1/2 have length £ in both a, r. Again, only one of the two triangulations can 
flip at either x ± 1/2. In the worst case scenario, x ± 1/2 can either both increase to £ + 2, 
or both decrease to £ — 2, or one can increase to £ + 2 and the other decrease to £ — 2. 
These produce respectively the terms 

U(r,r) = -%c,'^(l - a" 2 ) = i!!^ A(„,r); 

^^ = (r^ + a<rb)) A( "' T) ' 

Thus it suffices to show that, for any A < 1, we can find a > 1 and 5 > such that 

f 2aA 2 2 ] 

max \iTA2' oTiTa 2 ) ) sS1 ~' 5 - 

This is easily verified. □ 

Note: The bound in Theorem 15.31 should be compared with the known bound T m i x = 
0(n 3 logra) in the case of ID triangulations with A = 1 (and no constraints). This follows 
from a one-to-one correspondence between ID triangulations and lattice paths, and the 
results of Wilson [23] on the mixing time of the flip dynamics for such paths. 

6. Lower bounds on the mixing time for A > 1 

As a general rule, our lower bounds on the mixing time are obtained by exploiting the 
fact that it takes a long time to change the orientation of an initially long edge. Below, we 
say that an edge of a triangulation is positively oriented (respectively, negatively oriented) 
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if its leftmost endpoint lies below (respectively, above) its rightmost endpoint. A vertical 
or horizontal edge is neither positively nor negatively oriented. 

Our first result is an exponential lower bound on the mixing time (in the absence of 
constraints) that holds for any A > 1. This proves Theorem 11.41 in the Introduction. 

Theorem 6.1. For every A > 1, there exists c > such that the mixing time of the Glauber 
dynamics on triangulations o/A^ n in the absence of constraints satisfies T m ; x ^ exp (c (m + n)) 
for all m,n ^ 1. 

Proof For any A C £1 we write 

Z(A) = J>H 

for the partition function restricted to A. If dA denotes the set of a E A that are adjacent 
to A c in the flip graph, then the standard conductance bound (see, e.g., [TTJ Theorem 7.3]) 
will allow us to deduce the theorem once we find a set A such that n(A) ^1/2 and 

fg^e-^). (6.18) 

Assume w.l.o.g. that the horizontal coordinate n ^ m (else the same argument applies 
with the roles of n and m reversed). We define A C as the set of triangulations a £ 0, 
such that every internal midpoint with half-integer vertical coordinate v is not positively 
(resp., not negatively) oriented if v + ^ is odd (resp., even). Notice that any a G A has 
the "herringbone" structure illustrated in Figure [T2j In particular, any triangulation in A 
consists of ID triangulations in each horizontal layer. Moreover, if a E A then all edges 
with integer vertical coordinate are horizontal and frozen (unflippable). 
Now, the total length of any a £ A can be written as 

m 

\a\ = (m + l)n + (2n + l)m + V] L(a(j)), 

i=i 

where (m + l)n is the total length of all horizontal (unit length) edges, (2n + l)m is the 
sum of all vertical lengths of the edges of a, and L(a(j)) denotes the total horizontal 
length of the internal edges of the ID triangulation cr(j) in the j-th layer. 




FIGURE 12. Example of a triangulation of a 4 x 10 region, with the "her- 
ringbone" structure of the set A introduced in the proof of Theorem \6.1\ 



Below, we let Qf n denote the set of all ID triangulations such that all internal edges 
are not negatively oriented. From the well known map between lattice paths and ID 
triangulations, see, e.g., [8], one has that Qf n is in one-to-one correspondence with lattice 
paths (integer valued sequences with ±1 increments) of length 2n, which start at zero, 
end at zero, and stay non-negative. Moreover, under this correspondence, the total length 
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L(n) of 7T G is precisely the area under the path tt. The maximal element in Qf n , 
denoted tt, has total length 

n n— 1 

L(tt) = jN + J]i = n 2. 

1=1 1=1 

Clearly, er G (M implies that there exists a layer jf such that a(j) has an internal vertical 
edge, i.e., an edge with horizontal length zero. In particular, on that layer we must have 
L(a(j)) ^ n 2 — n. Therefore, 

E\LM 
,:<>,•, :/..: -)< ra 2 r, A 

< m J] A LW " n2 P„(A:)A- fc , 

iren+ n :L(n) < n 2 -n k^n 

where we have estimated T, 7T 'en+ > = A ™ 2 > and 

we use p n (k) to denote the 

number of tt G 0^ which have ^(vr) = n 2 — k. From the correspondence with lattice 
paths it is seen that p n (k) is bounded by the number p(k) of all partitions of the integer k. 
Using the well known fact that p(k) = e°(^^ we arrive at the desired bound (I6.18p . □ 

The lower bound in Theorem 16.11 can be improved as follows in the case m <C \fri. 

Theorem 6.2. For every A > 1, there exists c > such that the mixing time of the Glauber 
dynamics on triangulations o/A^ n , in the absence of constraints, satisfies T m j x ^ exp (cn 2 /■::: 
for all n ^ m 2 /c, m ^ 1. 

Proof. Take e G (0, 1) and divide the midpoints in A := A mjn into three regions as follows: 
Af is the set of x G A with horizontal coordinate between and en, A c is the set of 
x G A with horizontal coordinate between en and (1 — e)n, and A r is the remainder of A. 
Consider the set Ai C A c of all internal midpoints x G A c whose vertical coordinate is 
half- integer. Let A C f2 denote the set of triangulations a such that every edge a x with 
x G Ai with vertical coordinate v is positively or negatively oriented according to the 
parity of v + ^- I n particular, adjacent layers are required to have one-dimensional edges 
of opposite orientation throughout the whole central region. This is a relaxed version of 
the set A appearing in the proof of Theorem 16.11 (see Figure [T2l) . where the herringbone 
structure is now imposed only in the central region A c . 

We proceed by estimating Z(dA)/Z(A) as in the proof of Theorem 16. 11 Clearly, we can 
bound Z(A) from below by considering only the maximal one-dimensional triangulations 
in each layer j throughout the whole region A, and taking every edge with integer 
vertical coordinate to be horizontal (unit length), giving 

Z(A) ^ A mn2+ ° {mn) . 

Now let a G OA. The length of an edge whose midpoint is in A^ U A r is at most en + m. 
Since there are at most 8enm such midpoints, the total contribution of these edges to \o~\ 
is at most 8(en) 2 m + 8enm 2 . Moreover, there must be a midpoint xq G Ai whose edge o~ Xo 
is vertical (unit length) . The contribution of the layer containing xq to the total length of 
a is at most (un) 2 + ((1 — n)n) 2 , where e ^ u ^ 1 — e. All other layers contribute total 
length at most (m — l)n 2 . Hence for any a G dA, we have 

\a\ < 8(en) 2 m + 8enm 2 + (m — l)n 2 + (ii*n) 2 + ((1 — ti*)n) 2 , 
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where u* G [e, (1 — e)] maximizes (un) 2 + ((1 — u)n) 2 . Clearly, we have u* ~ e (or by 
symmetry u* ~ (1 — e)). Since there are at most 2 3mn triangulations in total [T], we get 



Z(dA) 
~Z(Aj 



yS(en) 2 m+d,enrn 2 -n 2 + {en) 2 +((l-e)n) 2 +0{mn) 



Thus taking e = c/m and n ^ m 2 /c for some sufficiently small constant c > (independent 
of n, m) gives 

g(g^) / -c'n 2 /m 

Z(A) - C 

This completes the proof. □ 



Note: While the above bound is tight in the case m = O(l), it becomes progressively 
weaker as m increases with n, and vacuous when m ~ y/n. In particular, it is far from 
the conjectured behavior T m ; x = e ^( mn ( m + n )) stated in Conjecture ll.il The deterioration 
of the bound as m increases is apparently an artifact of our proof, which is based on an 
essentially one-dimensional "herringbone" structure similar to that used in the proof of 
Theorem 16.11 

We end this section with a much weaker lower bound that holds for all A > 0. This 
shows in particular that the upper bound of Theorem 15.11 for the small A regime is tight. 

Proposition 6.3. There exists c > such that the mixing time of the Glauber dynamics 
on triangulations of A^ n in the absence of constraints satisfies T mix cmn(m + n) for 
all m, n 1 and all A > 0. 

Proof. Assume that n m. Let x = (1/2, n/2) £ A and let the initial condition a be 
an arbitrary triangulation in which a x is the edge from (0,0) to (l,n). Let A be the 
set of triangulations in which a x is not positively oriented. By symmetry, the stationary 
probability of A is at least 1/2 (with strict inequality if x is of type 1). Then, the total 
variation distance at time t can be bounded below as 

I|P*(<V)-HI >^-p\*,A) >^-F a (T A ^t), 

where ta denotes the hitting time of the set A, and P CT is the probability on trajectories of 
the Markov chain with initial condition a. Let A' C A be the set of midpoints (1/2, i) for 
i = 0, 1/2, 1, . . . , n. In order to flip o x all the way to an edge that is not positively oriented, 
we need to perform at least f2(n 2 ) flips of the edges with midpoints in A', as can be seen 
easily from the correspondence with lattice paths recalled in the proof of Theorem 16.11 
Thus P (T (' r A ^ t) is bounded above by the probability that out of t updates in A, at 
least Q(n 2 ) of them are in A'. Each update independently falls in A' with probability 

jj^- = 0[l/m). By Markov's inequality it follows that, for some constant C > 0, 



Taking t = mn 2 /AC yields T m ; x = fi(ran(m + nj). 



□ 
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7. Future Work 

Our results suggest a number of immediate open questions. 

(1) Can the spatial mixing property (Theorem \1.2\i and the polynomial mixing time 
bound (Theorem I1.3P be extended to the entire regime A < 1? This would in 
particular complete the verification of part (a) of Coniecture ll.il These questions, 
which are related, are likely to require a deeper investigation of the geometry of 
triangulations along the lines begun in Section SJ 

(2) Can the stronger exponential lower bound on mixing time in part (b) of Con- 
jecture 11.11 be proved in the regime A > 1? This will require the use of more 
sophisticated rigid triangulations than the "herringbone" structures used in Sec- 
tion [6l 

(3) The unweighted case A = 1 seems particularly challenging, as our results strongly 
suggest that it corresponds to a "phase transition." We conjecture that the spatial 
and temporal mixing properties here are likely to be very delicate. 

(4) In the "super-critical" regime A > 1, the model appears to exhibit various "phases" 
according to the direction of alignment of the edges. It would be interesting to 
describe these phases and their contributions to the equilibrium distribution. 
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